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ABSTRACT 

We describe observations of 9.7 /xm silicate features in 97 AGNs, exhibiting a 
wide range of AGN types and of X-ray extinction toward the central nuclei. We 
find that the strength of the silicate feature correlates with the HI column density 
estimated from fitting the X-ray data, such that low HI columns correspond to 
silicate emission while high columns correspond to silicate absorption. The be- 
havior is generally consistent with unification models where the large diversity in 
AGN properties is caused by viewing-angle-dependent obscuration of the nucleus. 
Radio-loud AGNs and radio-quiet quasars follow roughly the correlation between 
HI columns and the strength of the silicate feature defined by Seyfert galaxies. 
The agreement among AGN types suggests a high-level unification with similar 
characteristics for the structure of the obscuring material. We demonstrate the 
implications for unification models qualitatively with a conceptual disk model. 
The model includes an inner accretion disk (< 0.1 pc in radius), a middle disk 
(0.1-10 pc in radius) with a dense diffuse component and with embedded denser 
clouds, and an outer clumpy disk (10-300 pc in radius). 

Subject headings: galaxies: active — galaxies: nuclei 

1. Introduction 

The arrangement of material around the central supermassive blackhole (SMBH) in 
active galactic nuclei (AGNs) is crucial in unification schemes, where the active nuclei clas- 
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sified as type 2 arise from the obscuration of broad optical emission lines in type 1 ob- 
jects ( lAntonuccil Il993l ). The structure of the circumnuclear material may be critical to 
the growth of the SMBH and may influence the feedback from the AGN and its effects 
on the formation and evolution of the host galaxies. The infrared (IR) spectral energy 
distribution (SED) arises from the heating of surrounding material, and its shape should 
indica te how the material is org anized. Ther e are at least three differ e nt possibilities: corn- 



pact (IPier fc Krolik 



19921. 119931) exten ded (IGranato k Danese Il994l: iGranato et all Il997r ) 
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Nenkova et al 
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20021) struc- 



and cloudy (iKrolik fc Begelmanl 1 19861 : iRowan-RobinsonI 
ture models all fit the data well (assuming an extra component of star formation in the 
compact model). The silicate emission and absorption features now being observed with 
Spitzer (e.g., Weedman et al. 2005) can be used to probe these possibilities further. 



In this paper, we present Spitzer Infrared Spectrograph (IRS; iHouck et al.l |2004| ) ob- 
servations of 9.7 /ini silicate features in 93 AGNs, including e xamples from th e literature 
feiebenmorgen et al.l[2005l : iHao et al.l[2005l : Isturm et allboosl: lOde et al.ll2006h . We com- 



plement them with groundbased measurements for 4 more (IRoche et al. 



19911 1 . By combining 



these observations with X-ray data that explore the line of sight toward the central engine, 
we can gain new insights into the unification scheme. 



2. Sample and Data Reduction 



Our sample (see on-line Table 1) includes a variety of types of AGN: radio- loud quasars, 
FR II radio galaxies, Seyfert 1 and Seyfert 2 galaxies, broad absorption-line quasars (BALQs), 
low-ionization nuclear emission-line region (LINER) galaxies, and non-BAL radio-quiet quasars 
including optically selected Palomar- Green (PG) quasars and IR quasars as selected by the 
Two-Micron All Sky Survey (2MASS) from lSmith et al.l (120021 ). Except for the objects from 
the literature, quasars of all types and FR II radio galaxies are the most luminous objects 
in their parent samples (which include upper limits in redshift in their definitions) while 
Seyfert 1 and Seyfert 2 galaxies are deri v ed to have high brightness and a broad range in HI 
column density from [Turner fc Pounds! (119891 ) and iRisaliti et al.l ( 119991 ). respectively. The 



whole sample spans a range of HI column density from 10 cm up to 10 



cm ^. 



We focus 



on the sample of 85 objects with available HI column densities in this study, although all 
97 objects have available silicate data. X-ray spectra were retrieved from the Chandra data 
archive for 8 sources. The data red uction and the m easurement of column densities through 
power-law fits are as described in IShi et al.l (120051 ). The column dens i ties o f most of the 
2MASS QSOs are based on the X-ray hardness ratio from lWilkes et al. J2002h assuming an 



intrinsic power law X-ray spectrum with a photon index of 1.7. The associated uncertainty 
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is estimated as a factor of 3 corresponding to a change of ~1 in the photon index. The 
column density for the remaining sources is obtained from the hterature as shown in the 
on-hne Table 1. 

The IRS spectra presented for the first time in this paper, except for those of the Seyfert 
galaxies, were obtained using the standard staring mode. The intermediate products of the 
Spitzer Science Center (SSC ) pipeline Sll.0.2 were processed within the SMART software 
package (IHigdon et al.ll2004l ). The background was subtracted using associated spectra from 
the two nodded, off-source positions. This also subtracts stray light contamination from the 
peak-up apertures, and adjusts pixels with anomalous dark current. Pixels flagged by the 
SSC pipeline as "bad" were replaced with a value interpolated from an 8-pixel perimeter sur- 
rounding the suspect pixel. The spectra were extracted using a 6.0 pixel fixed-width aperture 
for the short low module (SL), and 5.0 pixels for the long low module (LL). The spectra were 
calibrated using a spectral response function derived from IRS spectra and Kurucz stellar 
models for a set of 16 Sun-like stars that exhibit: 1) high signal-to-noise observations, 2) no 
residual instrumental artifacts, and 3) no signs of IR excess. The absolute flux density scale 
is tied to calibrator stars observed by the IRS instrument team and r eferenced to calibrated 
stellar models provided by the SSC (see also iBouwman et al. II2006I ). The uncertainties in 
the final calibration are dominated by random noise. The relative flux calibra tion across the 
spectrum is accurate to ~ 1 — 2% (IHines et al. II2006I : iBouwman et al. II2006I ). 



The data reduction was slightly different for the Seyfert galaxies first presented in this 
paper. The whole Seyfert sample consists of a mixture of point-like and extended sources. 
They are at relatively low redshift {z ~0.01) and are lower luminosity AGNs, so care was 
taken to minimize the contribution of the host galaxies. For the extended sources, a fixed- 
width column extraction that is narrow enough to exclude most of the extended source 
would be most appropriate. However, narrow fixed-width extraction windows introduced 
artifacts into the extracted spectra. Instead, a narrow expanding-aperture extraction was 
performed on all sources, point-like or extended. This was deemed to be an acceptable 
compromise between rejecting the extended source component and introducing artifacts into 
the spectrum. The spectra were extracted using an expanding extraction aperture defined 
to be 4 pixels wide at the wavelengths of 6 nm and 12 urn for the second and the first 
order of SL, respectively, and of 16 fim and 27 /xm for the second and the first order of LL, 
respectively. At the wavelength of 10 /im, the extra ction aperture i s arou nd 3.3 pixels. A 
detailed description of the data reduction is given by iGorjian et al. I (120061 ). 



The silicate feature strengths were estimated as follows. Narrow emission lines were 
removed from the spectra and they were then smoothed to a resolution of 0.1 /im. The 
continuum was defined by using a spline interpolation between the blue and red ends of the 
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IRS spectral range. We defined (all in rest-frame wavelengths) the blue end as 5-7.5 /im, 
while the red end was defined as 13-14 fim for the 13 objects observed only with the SL 
module. For six BALQs at higher redshifts, the red end was chosen to be 10.5-13 /im. For 
the remainder of the sample, the red spectral end was defined as 25-30 /xm. For nine objects 
with the red end contaminated by silicate emission (Figure [T](b)), the continuum at this end 
was estimated to be below that observed by <20% at the reddest point, based on the spectra 
of objects without contamination. For twelve objects with the blue end contaminated by 
the strong 6-8 /im aromatic feature (Figure [D^c)), the continuum in this bandpass was given 
by two segments between 5 and 5.5 fim where the contamination is negligible. The fit was 
judged to be good when it matched the flux at both ends within the noise and the curvature 
of the continuum varied gradually from one end to the other over the silicate region. The 
uncertainty in continuum fitting was obtained by adjusting the flux of the continuum up 
and down over the fitting spectral range. For the objects with continuum not contaminated 
by either aromatic or silicate features as shown in Figure [T](a), we adjusted the fit until the 
continuum was just above (for upper error) or below (for lower error) the observed flux over 
the whole spectral fitting range. To estimate the uncertainty for the objects with red end 
contaminated by the silicate feature (Figure [T](b)), we first fixed the blue end and adjusted 
the red end continuum up to the red end point that matched the upper-envelope of the 
underlying continuum plus noise and down by the same amount as an estimate of the errors 
of the continuum fit. To account for the effects of noise in the blue end, we added another 4% 
uncertainty (the mean value for objects without aromatic or silicate feature contamination). 
For the objects with blue end contaminated by the aromatic feature (Figure [T](c)), we fixed 
the red end and adjusted the blue end to get the uncertainty in a similar way where the 
upper-limit of the underlying continuum is the flux at 7 fim (the median wavelength between 
the peaks of the 6.2 /im and 7.7 /im aromatic features). Another 4% uncertainty was added 
to this uncertainty to account for the effects of noise in the red end. There are no objects 
with red end contaminated by the silicate feature and blue end contaminated by the aromatic 
features. 

The strength of the 9.7 /im silicate feature is defined as (Ff — Fc)/Fc, where F{ and Fc are 
the observed flux density and underlying continuum flux density, respectively, at the peak 
(for emission) or the minimum (for absorption) of the silicate feature. The feature strength 
in this definition is a direct measure of the optical depth for the silicate absorption. Figure [2] 
shows the IRS spectra presented for the first time in this study in order of HI column density. 
No ice or hydrocarbon absorption is found in Figure [2l Especially at 5-8 /im, where there is 
no contamination of the silicate feature, the spectra are well described by power laws except 
for those with aromatic emission. Therefore, once the continuum had been fitted, we could 
calculate the feature strengths unambiguously. 
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3. Results 

Figure [3] shows the strength of the 9.7 /im sihcate feature as a function of HI column den- 
sity. The right y-axis of Figure [3] shows the IR-absorbing column density es timated from the 



silicate absorption feature by assumin g T<^,7j,,m. = Inf-Fr/Ff), 74v Mq 7^,^=19 (IRoche fc Aitken 



19851 ) and A^/A^j^=0.62xl0-2icm-2 flSavage fc MathisI Il979r ). The silicate feature varies 



from emission (+) to absorption (-) as the X-ray spectra become more heavily obscured. 
The trend is also demonstrated by the composite spectra in different bins of HI column as 
shown in Figure HI Since the silicate feature is broad, the co mposite spectra are the ge omet- 



ric mean spectra to conserve the global spectral shape (e.g. IVanden Berk et al.ll200ll ). The 
relationship between HI column and silicate feature behavior is generally consistent with the 
AGN unification scheme where material surrounding the central SMBH obscures both the 
X-ray and silicate emissions. 

The correlation (black lines in Figure [3]) is defined by all types of AGN. To investigate 
the behavior for individual AGN types, we unified Seyfert 1 and Seyfert 2 galaxies as one type 
of AGN (Seyfert galaxies). Similarly, FR II radio galaxies and radio- loud QSOs are classified 
as radio-loud AGNs while PG and 2MASS QSOs are classified as radio-quiet QSOs. Table 2 
shows the linear fits to these AGN types where the limits to the HI column density are treated 
as detections during the fitting. The fits for different AGN types are characterised by large 
uncertainty and are generally consistent with each other within the uncertainty. Radio-quiet 
QSOs (2MASS and PG QSOs) exhibit a small (less than two standard deviations) deviation 
from other two types of AGN. As shown in Figure [3], this deviation arises because the PG 
QSOs with low HI columns have larger silicate strengths. 

To quantify the discrepancy between different types of AGN, we examined the proba- 
bility that the other AGN types follow the Seyfert correlation by producing a Monte Carlo 
theoretical distribution with 10000 data points using the Seyfert correlation (red line in Fig- 
ure [3]) with associated scatter. Table 3 shows the result of three tests: 1.) the probability 
from a K-S test that the given AGN type has the same distribution as the theoretical distri- 
bution; 2.) the significance that two distributions have the same mean (indicated by a value 
greater than 0.05); and 3.) the significance that two distributions have the same variance 
(indicated by a value greater than 0.05). The high K-S probability for the Seyfert galaxies 
indicates that the theoretical distribution is a good representative of them (as it should be). 
As shown in Table 3 and Figure [31 radio-loud AGNs (radio-loud QSOs and FR II galaxies) 
are consistent with the Seyfert correlation and radio-quiet QSOs (2MASS and PG QSOs) 
show slightly higher silicate strengths at low HI column compared to Seyfert galaxies. 



Sturm et al.l (120061 ) recently found that six type 2 QSOs in the HI column range of 10 



21.5 



10 cm do not show any silicate feature, slightly higher but still consistent with those of 
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the Seyfert galaxies at the same range of HI columns. As shown in Figure [3l the three BALQs 
and the one LINER in this study deviate from the Seyfert correlation significantly. More 
objects are needed to address whether these two types follow the Seyfert correlation. Given 
that most PG QSOs have upperlimit measurements of HI columns, the intrinsic deviation of 
radio-quiet QSOs from the Seyfert trend should be smaller than we calculate. Although the 
intrinsic slope of the Seyfert correlation should be smaller due to the lowerlimit measurement 
of the HI column for half the Seyfert 2 galaxies, the effect on the comparision is small because 
it is dominated by the column range where Seyfert galaxies have detected HI columns. In 
Table 2, we also list the linear fits to the Seyferts and the whole sample excluding objects 
with limit measurements of the HI columns. The comparision between the new slopes shows 
that the consistency between different AGN types may become even better. 

Therefore, except for LINERs and BALQs, the remaining non-Seyfert AGN types follow 
more or less the Seyfert correlation. The rough agreement among AGN types indicates 
that the geometry of the circumnuclear material and our viewing angle are the primary 
factors influencing the relation between the silicate feature and the X-ray attenuation. This 
correlation is nonetheless curious because, as shown in Figure [3l the column required for the 
observed X-ray obscuration levels is up to ~10^^ cm~^, two orders of magnitude larger than 
that required to produce the silicate absorption, ~ 10^^ cm~^. This implies that along many 
lines of sight there must be X-ray absorbing material that is not contributing to the silicate 
absorption. 

The dispersion in the correlation is much larger than the uncertainties in silicate-feature 
strength and HI column density. The large dispersion appears to be characteristic of all the 
AGN types with large numbers of objects (Seyfert galaxies, radio-quiet QSOs and radio- 
loud AGNs). This is consistent with a common mechanism (the circumnuclear geometry) 
regulating the correlation and dispersion in Figure [31 Due to different physical aperture 
diameters from 500 pc to galaxy-scale in size, the differing amounts of extended emission 
from star-forming regions may contribute to the scatter. However, such contributions must 
be small. As shown in Figure [2l twelve of 85 objects have strong aromatic features and their 
silicate features may be contaminated by star-forming regions. However, the dispersion of the 
correlation is almost the same for the 73 objects without aromatic features. Another possible 
contribution to the scatter is the gas-to-dust ratio. However, as indicated by Figure [31 the 
gas-to-dust ratio needs to vary nearly three orders of magnitude to account for the scatter 
in the correlation. Such large variation is unreasonable given the similar IR SEDs and IR 
luminosities for objects with similar HI column but different silicate features as shown in 
Figured for example, NGC 4941 versus NGC 328L 

Another characteristic of the correlation is that several Compton-thick AGNs have sil- 
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icate absorptions that are much weaker than predicted. Again, this imphes that additional 
absorbing material must be present in the Compton-thick sources, but placed so it does not 
obscure the IR emission. 



Discussion 



We now describe a conceptual model to explain the large difference in the absorbing 
columns for the X-rays and infrared, but that can also account for the correlation in Figure [31 
The goal is to find a geometry for the circumnuclear material that can explain Figure [3] solely 
in terms of variations in viewing angle. As mentioned in the Introduction, there are three 
classes of model that ar e generally successfu l in fitting the IR SEDs of specific AGN types: 
1.) a compa ct disk (e.g. iPier fc Kroliklll992l ): 2.) a more extended disk (radius of hundreds 



of pc) (e.g. iGranato et al.l 119971 ): and 3.) a disk with clumps or clouds of material (e.g. 



Nenkova et al.ll2002l ). We now explore whether the ideas behind these three types of model 
can be combined to provide a possible explanation for the behavior of the silicate feature. It 
is beyond the scope of this paper to compute a quantitative model to fit the data. However, 
by confining our argument to combinations of features in models already shown to fit other 
aspects of AGN behavior, it is likely that our explanation of the silicate-X-ray correlation 
will also be compatible with the other observations of AGNs. Overall, our model r equires 

2002h and a 



Nenkova et a. 



Pier fc Krolik 



both a component of material similar to the cloudy model (e.g. 
diffus e component with outer radius similar to the compact disk-model (e.g. 
1992) and with column density similar to the model of the extended disk (e.g. iGranato et al. 
1997( 1. 



In Figure O we show the hypothetical disk geometry. We include an inner accretion disk 
(AD), and in the same plane a middle disk (MD) with a diffuse component (grey) joining 
to the AD. Denser clouds or clumps are embedded in this MD. It merges with a cloudy or 
clumpy outer disk (OD). The X-ray and UV radiation from the AD heat the dust in the MD 
and OD to produce IR emission and ionize the broad-emission-line clouds and the narrow- 
emission-line (NEL in Figure [5]) clouds, which are not shielded by the diffuse component. 
Compton-thick X-ray obscuration can arise either in the AD, or in the clouds in the MD. 
Irregularities in the AD or the passage of clouds through the line of sight may be responsible 



for the variations i n X-ray obscuration observed toward some AGNs (e.g. lElvis et al.l 12004 



Risaliti et al.ll2005l ). However, for a model with only the central disk and the clouds, there 
is a large probability of viewing the central X-ray emission directly without any extinction 
even for an inclined disk (along the equatorial plane). This is inconsistent with Figure [3l 
where the lower-envelope of the data distribution shows a strong trend that the HI column 
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density increases with the depth of the sihcate absorption. The diffuse component of the 
MD is required to obscure the X-ray emission even if the hne of sight does not cross any 
cloud. This component is also the source of the silicate emission. The silicate absorption 
results from the material in the OD. 

We argue that the outer edge of the MD should not extend beyo nd around 10 pc and 



the dust within it can have sufficiently high temperature (> 300 K; iLaor fc Draindll993l ) 
to produce the silicate emission. Given the maximum silicate absorption of lO^^cm"^ and 
the maximum HI column of >10^^cm~^, the clouds in the MD should have A'^^^ > 10^^ cm~^ 
and produce significant obscuration of the X-ray emission for an intercepting line of sight 
but not obscure the silicate emission significantly due to a small covering factor. For a given 
strength of the silicate feature along a line of sight, the large variation in the HI column 
is caused by the variation in the number or column density of clouds that the line of sight 
intercepts and the minimum X-ray obscuration is due only to the diffuse component. Based 
on this concept, as the line of sight varies such that the silicate strength ranges from 0.0 
to -0.8 (the minimum value in the sample) as shown in Figure [3l the column density of 
the diffuse component varies from 10^^ cm~^ to 10^^ c m~^. The evidenc e for the cloudy 



QD is that the mid-IR images of NGC 1068 at 10 /xm flJaffe et all l2004f ) and NGC 4151 



(iRadomski et al.ll2003l ) at ~10 and ~18 /im show that the radius of any diffuse component 
should be smaller than 2 pc and 35 pc, respectively. However, the dust in the OD is relatively 
cool and emits any reprocessed energy mainly at far-IR wavelengths and thus may be missed 
for observations at mid-IR wavelengths. Since both the MD and OD are unresolved in the 
IRS beam, the observed silicate absorption can be provided by several clouds in the OD and 
the depth of the absorption feature is determined by the average number of clouds along a 
line of sight. This is because a cloud in the OD is not large enough to cover the whole MD. 



5. Conclusions 

We report observations of 9.7/im silicate features in 97 AGNs. The features vary from 
emission to absorption with increasing HI column density, consistent with unification mod- 
els. Radio-loud AGN (radio-loud QSOs and FR II galaxies) and radio-quiet QSOs (PG and 
2MASS QSOs) lie roughly on the Seyfert-correlation between HI column and silicate feature 
strength. The behaviors of LINERs and BALQs are not clear due to the small number of 
objects for these two AGN types. The scatter in the relation is large and several Compton- 
thick AGNs do not show deep silicate absorption. Qualitatively, the correlation requires a 
circumnuclear disk geometry with an accretion disk outside of which is a middle disk with 
high density and with even denser clouds embedded (0.1-10 pc in radius), co-ahgned with 
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an outer clumpy disk (10-300 pc in radius). The similarity of the behavior of various types 
of AGN suggests that this disk geometry may be typical for AGN in general. 

We thank the anonymous referee for detailed comments. We also thank Roberto 
Maiohno for helpful suggestions. Support for this work was provided by NASA through 
contract 1255094 issued by JPL/ Cahfornia Institute of Technology. 
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Table 1. Source Characteristics 



Sources 


z 


Type 


Strength 


Reference 


N§ 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



PG0050+124 0.058 PG 0.38 2 

PG0052+251 0.155 PG O.SSto.oe 

PG0804+761 0.100 PG 0.60 2 

PG0953+414 0.234 PG 0.40to.o8 



Note. — Column (1): The sources. Column (2): Redshift. Column (3): The types 
of AGNs. 'PG': PG quasar; 'RLQ': radio-loud quasar; 'Syl': Seyfert 1 galaxies; 
'Sy2': Seyfert 2 galaxies; 'BALQ': broad absorption-line quasar; 'FRII': FR II ra- 
dio galaxies; 'LINER': low-ionization nuclear emission-line region; '2MQ': 2MASS 
quasar. Column(4): The strength of the silicate feature a s defined in § 2. Column (5) : 
Th e reference s for t he si licate data: (0) Th is wo rk; fl) ISiebe nmorge n et al. | (| 2005|): 
(2) iHao et al.l ^m^; f3') ISturm et al.l ll2005h ; (41 iRoche et al.l l ll99j) ; (5) lOgle et al I 



Column (6): The intrinsic HI column density in the unit of 10^ 



Col- 



COOi 

umn (7): References for t h e HI colum n densities: (6 ) This work; (7) Donato^^al 
1 20031 ); (8) iBrunner et al.l lll997l') : (9) IPorguet et al.l (I20q4l'l: ( 10) Reeves fc Turner 




Leighlv et aU 11997 
Dewangai^_etal. j2003l) ; 
Wilkes et al.lll2002h; (18 



Gallagher et al.l lll999h " (21 

" ' 200d) : (24) lls"obe et al. 



151 



Gallagher et al.l (|2002|): ( 13) iMathur et al 



Turner &: Poundi 
19961) 



I lll989l) 



(161 



Risaliti et al 



Wang et al. ( 
Wilkes et al.l 1I2OO5I ) 



Belsole et al.l 



ll2002h 



(22) 



19)|Ptak et al.lll2004) ; (20) 
Worrall et 3111120011 ) ;(23) 



The complete version of this table is in the electronic edition of the Journal. 
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Table 2. The linear fits to objects with 
different AGN types 



Type 


A 


B 


Number 


(1) 


(2) 


(3) 


(4) 


Seyfert 


2.6±0.7 


-0.12 ± 0.03 


38 


2MQ & PG 


5.5±1.2 


-0.25 ± 0.06 


28 


FRII & RLQ 


3.9±1.6 


-0.17 ± 0.07 


15 


all objects 


3.3±0.5 


-0.15 ± 0.02 


85 


Seyferti 


4.6±1.3 


-0.21 ± 0.06 


23 


all objects^ 


5.1±0.9 


-0.23 ± 0.04 


61 



Note. — Column(l): The AGN types. See caption to Figure 
3. Column(2) and Column(3): Parameters of the linear fit S = 
A + B * log{N^) where S is the strength of the silicate feature 
as defined in § 2 and A''^ is the HI column density in the unit of 
cm~^. Both parameters A and B are unitless due to the definition 
of the silicate strength. Column(4): the total number of objects 
for each fit. 

^Fits excluding objects without measured HI columns (i.e., upper 
and lower limits to the HI are excluded). 
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Table 3. Tests of the correlation 



Type 


Num 


K-S Prob.{%) 


Mean 


Variance 


(1) 


(2) 


(3) 


(4) 


5 


Sy 


38 


90 


1.0 


0.85 


2MQ & PG 


28 


2 


0.02 


0.35 


RLQ & FR II 


15 


10 


0.06 


0.40 



Note. — Column{l): The AGN type. See caption to Figure [3] 
Column(2): the total number of a given AGN type. Column(3): 
The probabihty from the K-S test that the given AGN type has 
the same distribution as a theoretical Monte Carlo distribution 
produced by the correlation defined by the Seyfert galaxies. Col- 
umn (4) : The significance that the distribution of the given AGN 
type and the theoretical distribution have the same mean. A value 
greater than 0.05 indicates the same mean for the two distribu- 
tions. Column (5): The significance that the distribution of the 
given AGN type and the theoretical distribution have the same 
variance. A value greater than 0.05 indicates the same variance 
for the two distributions. 
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Fig. 1. — IRS spectra of 3 objects with fits. Tlie daslied line is the fitted continuum and 
the dotted hne is the estimated uncertainty in the continuum fitting. The scahng factors for 
(a) MCG-2-58-22, (b) Fairall9 and (c) NGC4388, are 30, 5 and 1.5, respectively. 
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Fig. 2. — The IRS spectra for objects with available X-ray data in order of the HI column 
density. Narrow emission lines have been removed from the spectra (see text). The dashed 
lines are the fitted continua and the dotted lines show the estimated uncertainties in the 
continuum fitting. The scaling factor to normalize the spectra at 5 /^m and the logarithm of 
the HI column density in units of cm~^ are given for each object. 
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Fig. 2. — Continued 
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Fig. 2. — Continued 
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Fig. 2. — Continued 




Fig. 2. — Continued 



Rest-frome wavelength (yu-m) 



- 21 - 



1.5[ ' 




-1.0 [ \ , \ , \ , \ , ] 

1018 1020 1022 1024 1026 

Nh [cm-2] 

Fig. 3. — The strength of the sihcate feature as a function of HI column density. The 
strength of the sihcate feature is defined as (Ff — Fc)/Fc, where Ff and Fc are the observed 
flux density and underlying continuum flux density, respectively, at the peak (for emission) 
or the minimum (for absorption) of the silicate feature. The black line is the linear fit to all 
objects while the red hne is the fit to Seyfert galaxies. 'PC: PG quasar; 'RLQ': radio-loud 
quasar; 'Syl': Seyfert 1 galaxies; 'Sy2': Seyfert 2 galaxies; 'BALQ': broad absorption-line 
quasar; 'FRII': FR II radio galaxies; 'LINER': low-ionization nuclear emission-line region; 
'2MQ': 2MASS Quasar. 




Fig. 4. — The composite spectra of AGNs in different HI column bins. The number in 
parenthesis is the number of objects used for the composite spectrum in each bin. 



-23- 




^ ^ ^ ^ ^ ^ 

0.1 1 10 100 

Equatorial plane[(L45) V2 pc] 



Fig. 5. — The structure of the material surrounding the central blackhole in the first quarter 
section. The whole structure is symmetric about the disk axis and the equatorial plane. Prom 
inside to outside: 1.) the inner accretion disk (AD), which produces X-ray and UV radiation 
ionizing the narrow-emission-line (NEL) and broad-emission-line clouds, and heating the 
dust; 2.) the middle disk (MD) with a diffuse component (grey) and with denser embedded 
clouds - the diffuse component produces the silicate emission while the embedded clouds 
heavily obscure the central X-ray emission when the line of sight intercepts them; and 3.) the 
outer disk (OD) with clouds that obscure the silicate and X-ray emission, and are responsible 
for the far-IR emission. Four fines of sight indicate: 1.) sificate emission with low HI column; 
2.) silicate emission with high HI column; 3.) silicate absorption with low HI column; 4.) 
silicate absorption with high HI column. 



